Serial adaptive laboratory evolution enhances mixed carbon metabolic capacity of Escherichia coli.

Microbes have inherent capacities for utilizing various carbon sources, however they often exhibit sub-par fitness due to low metabolic efficiency. To test whether a bacterial strain can optimally utilize multiple carbon sources, Escherichia coli was serially evolved in L-lactate and glycerol. This yielded two end-point strains that evolved first in L-lactate then in glycerol, and vice versa. The end-point strains displayed a universal growth advantage on single and a mixture of adaptive carbon sources, enabled by a concerted action of carbon source-specialists and generalist mutants. The combination of just four variants of glpK, ppsA, ydcI, and rph-pyrE, accounted for more than 80% of end-point strain fitness. In addition, machine learning analysis revealed a coordinated activity of transcriptional regulators imparting condition-specific regulation of gene expression. The effectiveness of the serial adaptive laboratory evolution (ALE) scheme in bioproduction applications was assessed under single and mixed-carbon culture conditions, in which serial ALE strain exhibited superior productivity of acetoin compared to ancestral strains. Together, systems-level analysis elucidated the molecular basis of serial evolution, which hold potential utility in bioproduction applications.

Reduction-to-synthesis: the dominant approach to genome-scale synthetic biology.

Advances in systems and synthetic biology have propelled the construction of reduced bacterial genomes. Genome reduction was initially focused on exploring properties of minimal genomes, but more recently it has been deployed as an engineering strategy to enhance strain performance. This review provides the latest updates on reduced genomes, focusing on dual-track approaches of top-down reduction and bottom-up synthesis for their construction. Using cases from studies that are based on established industrial workhorse strains, we discuss the construction of a series of synthetic phenotypes that are candidates for biotechnological applications. Finally, we address the possible uses of reduced genomes for biotechnological applications and the needed future research directions that may ultimately lead to the total synthesis of rationally designed genomes.

Minireview: Engineering evolution to reconfigure phenotypic traits in microbes for biotechnological applications.

Adaptive laboratory evolution (ALE) has long been used as the tool of choice for microbial engineering applications, ranging from the production of commodity chemicals to the innovation of complex phenotypes. With the advent of systems and synthetic biology, the ALE experimental design has become increasingly sophisticated. For instance, implementation of in silico metabolic model reconstruction and advanced synthetic biology tools have facilitated the effective coupling of desired traits to adaptive phenotypes. Furthermore, various multi-omic tools now enable in-depth analysis of cellular states, providing a comprehensive understanding of the biology of even the most genomically perturbed systems. Emerging machine learning approaches would assist in streamlining the interpretation of massive and multiplexed datasets and promoting our understanding of complexity in biology. This review covers some of the representative case studies among the 700 independent ALE studies reported to date, outlining key ideas, principles, and important mechanisms underlying ALE designs in bioproduction and synthetic cell engineering, with evidence from literatures to aid comprehension.

Genome-wide CRISPRi screen identifies enhanced autolithotrophic phenotypes in acetogenic bacterium Eubacterium limosum.

Acetogenic bacteria are a unique biocatalyst that highly promises to develop the sustainable bioconversion of carbon oxides (e.g., CO and CO2) into multicarbon biochemicals. Genotype-phenotype relationships are important for engineering their metabolic capability to enhance their biocatalytic performance; however, systemic investigation on the fitness contribution of individual gene has been limited. Here, we report genome-scale CRISPR interference screening using 41,939 guide RNAs designed from the E. limosum genome, one of the model acetogenic species, where all genes were targeted for transcriptional suppression. We investigated the fitness contributions of 96% of the total genes identified, revealing the gene fitness and essentiality for heterotrophic and autotrophic metabolisms. Our data show that the Wood-Ljungdahl pathway, membrane regeneration, membrane protein biosynthesis, and butyrate synthesis are essential for autotrophic acetogenesis in E. limosum. Furthermore, we discovered genes that are repression targets that unbiasedly increased autotrophic growth rates fourfold and acetoin production 1.5-fold compared to the wild-type strain under CO2-H2 conditions. These results provide insight for understanding acetogenic metabolism and genome engineering in acetogenic bacteria.

Revealing Causes for False-Positive and False-Negative Calling of Gene Essentiality in Escherichia coli Using Transposon Insertion Sequencing.

The massive sequencing of transposon insertion mutant libraries (Tn-Seq) represents a commonly used method to determine essential genes in bacteria. Using a hypersaturated transposon mutant library consisting of 400,096 unique Tn insertions, 523 genes were classified as essential in Escherichia coli K-12 MG1655. This provided a useful genome-wide gene essentiality landscape for rapidly identifying 233 of 301 essential genes previously validated by a knockout study. However, there was a discrepancy in essential gene sets determined by conventional gene deletion methods and Tn-Seq, although different Tn-Seq studies reported different extents of discrepancy. We have elucidated two causes of this discrepancy. First, 68 essential genes not detected by Tn-Seq contain nonessential subgenic domains that are tolerant to transposon insertion, which leads to the false assignment of an essential gene as a nonessential or dispensable gene. These genes exhibited a high level of transposon insertion in their subgenic nonessential domains. In contrast, 290 genes were additionally categorized as essential by Tn-Seq, although their knockout mutants were available. The comparative analysis of Tn-Seq and high-resolution footprinting of nucleoid-associated proteins (NAPs) revealed that a protein-DNA interaction hinders transposon insertion. We identified 213 false-positive genes caused by NAP-genome interactions. These two limitations have to be considered when addressing essential bacterial genes using Tn-Seq. Furthermore, a comparative analysis of high-resolution Tn-Seq with other data sets is required for a more accurate determination of essential genes in bacteria. IMPORTANCE Transposon mutagenesis is an efficient way to explore gene essentiality of a bacterial genome. However, there was a discrepancy between the essential gene set determined by transposon mutagenesis and that determined using single-gene knockout strains. In this study, we generated a hypersaturated Escherichia coli transposon mutant library comprising approximately 400,000 different mutants. Determination of transposon insertion sites using next-generation sequencing provided a high-resolution essentiality landscape of the E. coli genome. We identified false negatives of essential gene discovery due to the permissive insertion of transposons in the C-terminal region. Comparisons between the transposon insertion landscape with binding profiles of DNA-binding proteins revealed interference of nucleoid-associated proteins to transposon insertion, generating false positives of essential gene discovery. Consideration of these findings is required to avoid the misinterpretation of transposon mutagenesis results.

Comprehensive 16S rRNA and metagenomic data from the gut microbiome of aging and rejuvenation mouse models.

The gut microbiota is associated with the health and longevity of the host. A few methods, such as fecal microbiota transplantation and oral administration of probiotics, have been applied to alter the gut microbiome and promote healthy aging. The changes in host microbiomes still remain poorly understood. Here, we characterized both the changes in gut microbial communities and their functional potential derived from colon samples in mouse models during aging. We achieved this through four procedures including co-housing, serum injection, parabiosis, and oral administration of Akkermansia muciniphila as probiotics using bacterial 16 S rRNA sequencing and shotgun metagenomic sequencing. The dataset comprised 16 S rRNA sequencing (36,249,200 paired-end reads, 107 sequencing data) and metagenomic sequencing data (307,194,369 paired-end reads, 109 sequencing data), characterizing the taxonomy of bacterial communities and their functional potential during aging and rejuvenation. The generated data expand the resources of the gut microbiome related to aging and rejuvenation and provide a useful dataset for research on developing therapeutic strategies to achieve healthy active aging.

Systems Biology on Acetogenic Bacteria for Utilizing C1 Feedstocks.

With a presence of the Wood-Ljungdahl pathway, acetogenic bacteria are capable of converting C1 feedstocks into biomass and various metabolites, receiving industrial interest in microbial production of biochemicals derived from C1 substrates. To understand C1 feedstock fermentation using acetogenic bacteria, most of the studies have focused on revealing their carbon assimilation and energy conservation systems. Despite the determination of the essential mechanisms, a fundamental understanding of acetogenic bacteria and the associated complex regulatory systems remains unclear and is needed for rational strain design. For this purpose, systems biology is a suitable approach for investigating genome, transcription, translation, regulation systems, and metabolic flux, providing a glimpse of the relationship between the genotype and phenotype of the organisms. This chapter will cover recent systems biology applications on acetogenic bacteria and discuss the cellular responses during C1 feedstock fermentation along with the regulatory systems that orchestrate cellular processes.

Functional role of a novel algicidal compound produced by Pseudoruegeria sp. M32A2M on the harmful algae Alexandrium catenella.

A marine phytoplankton dinoflagellate, Alexandrium sp. is known to cause worldwide harmful algal blooms, resulting in paralytic shellfish poisoning. In this study, we isolated a novel compound secreted by the marine bacterium Pseudoruegeria sp. M32A2M, and showed that it displays algicidal activity against A. catenella (group I). The molecular structure of the compound was analyzed by using 1H nuclear magnetic resonance (NMR), 13C NMR, and gas chromatography-mass spectrometry, which revealed that the compound was a diketopiperazine, cyclo[Ala-Gly]. Cyclo[Ala-Gly] induced a rapid decrease in the active chlorophyll a content and maximal quantum yield of photosystem II, leading to membrane disintegration after 24 h of its treatment. It showed the highest algicidal effect against diketopiperazines and also showed specific algicidal activities against several dinoflagellate species, but not for diatom species. In particular, cyclo[Ala-Gly] caused the transcriptional downregulation of the photosynthesis-related membrane complex in A. catenella, but not in the diatom Chaetoceros simplex. Based on structural modeling, we elucidated that cyclo[Ala-Gly] has a structure similar to that of plastoquinone, which transfers electrons by binding to the photosystem II core proteins PsbA and PsbD. This suggests a novel role for cyclo[Ala-Gly] as a potential inhibitor of photosynthesis.

Development of highly characterized genetic bioparts for efficient gene expression in CO2-fixing Eubacterium limosum.

Acetogenic bacteria demonstrate industrial potential for utilizing carbon dioxide (CO2) for biochemical production using the Wood-Ljungdahl pathway. However, the metabolic engineering of acetogenic bacteria has been hampered by the limited number of available genetic bioparts for gene expression. Here, we integrated RNA sequencing, ribosome profiling, differential RNA sequencing, and RNA 3'-end sequencing results of Eubacterium limosum to establish genetic bioparts, such as promoters, 5' untranslated regions, and transcript terminators, to regulate transcriptional and translational expression of genes composing of biosynthetic pathways. In addition, a transformation method for the strain was developed to efficiently deliver the obtained genetic bioparts into cells, resulting in a transformation efficiency of 2.5 × 105 CFU/µg DNA. Using this method, the genetic bioparts were efficiently introduced, and their strengths were measured, which were then applied to optimize the heterologous expression of acetolactate synthase and acetolactate decarboxylase for non-native biochemical acetoin production. The strategy developed in this study is the first report on integrating multi-omics data for biopart development of CO2 or syngas utilizing acetogenic bacteria, which lays a foundation for the efficient production of biochemicals from CO2 or syngas as a carbon feedstock under autotrophic growth conditions.

Synthetic 3'-UTR valves for optimal metabolic flux control in Escherichia coli.

As the design of genetic circuitry for synthetic biology becomes more sophisticated, diverse regulatory bioparts are required. Despite their importance, well-characterized 3'-untranslated region (3'-UTR) bioparts are limited. Thus, transcript 3'-ends require further investigation to understand the underlying regulatory role and applications of the 3'-UTR. Here, we revisited the use of Term-Seq in the Escherichia coli strain K-12 MG1655 to enhance our understanding of 3'-UTR regulatory functions and to provide a diverse collection of tunable 3'-UTR bioparts with a wide termination strength range. Comprehensive analysis of 1,629 transcript 3'-end positions revealed multiple 3'-termini classes generated through transcription termination and RNA processing. The examination of individual Rho-independent terminators revealed a reduction in downstream gene expression over a wide range, which led to the design of novel synthetic metabolic valves that control metabolic fluxes in branched pathways. These synthetic metabolic valves determine the optimal balance of heterologous pathways for maximum target biochemical productivity. The regulatory strategy using 3'-UTR bioparts is advantageous over promoter- or 5'-UTR-based transcriptional control as it modulates gene expression at transcription levels without trans-acting element requirements (e.g. transcription factors). Our results provide a foundational platform for 3'-UTR engineering in synthetic biology applications.

Genome-scale analysis of genetic regulatory elements in Streptomyces avermitilis MA-4680 using transcript boundary information.

BACKGROUND: The gram-positive bacterium, Streptomyces avermitilis, holds industrial importance as the producer of avermectin, a widely used anthelmintic agent, and a heterologous expression host of secondary metabolite-biosynthetic gene clusters. Despite its industrial importance, S. avermitilis' genome organization and regulation of gene expression remain poorly understood. In this study, four different types of Next-Generation Sequencing techniques, including dRNA-Seq, Term-Seq, RNA-Seq and ribosome profiling, were applied to S. avermitilis to determine transcription units of S. avermitilis at a genome-wide level and elucidate regulatory elements for transcriptional and translational control of individual transcription units. RESULT: By applying dRNA-Seq and Term-Seq to S. avermitilis MA-4680, a total of 2361 transcription start sites and 2017 transcript 3'-end positions were identified, respectively, leading to determination of 1601 transcription units encoded in S. avermitilis' genome. Cataloguing the transcription units and integrated analysis of multiple high-throughput data types revealed the presence of diverse regulatory elements for gene expression, such as promoters, 5'-UTRs, terminators, 3'-UTRs and riboswitches. The conserved promoter motifs were identified from 2361 transcription start sites as 5'-TANNNT and 5'-BTGACN for the - 10 and - 35 elements, respectively. The - 35 element and spacer lengths between - 10 and - 35 elements were critical for transcriptional regulation of functionally distinct genes, suggesting the involvement of unique sigma factors. In addition, regulatory sequences recognized by antibiotic regulatory proteins were identified from the transcription start site information. Analysis of the 3'-end of RNA transcript revealed that stem structure formation is a major determinant for transcription termination of most transcription units. CONCLUSIONS: The transcription unit architecture elucidated from the transcripts' boundary information provides insights for unique genetic regulatory mechanisms of S. avermitilis. Our findings will elevate S. avermitilis' potential as a production host for a diverse set of secondary metabolites.

Adaptive laboratory evolution of Escherichia coli W enhances gamma-aminobutyric acid production using glycerol as the carbon source.

The microbial conversion of glycerol into value-added commodity products has emerged as an attractive means to meet the demands of biosustainability. However, glycerol is a non-preferential carbon source for productive fermentation because of its low energy density. We employed evolutionary and metabolic engineering in tandem to construct an Escherichia coli strain with improved GABA production using glycerol as the feedstock carbon. Adaptive evolution of E. coli W under glycerol-limited conditions for 1300 generations harnessed an adapted strain with a metabolic system optimized for glycerol utilization. Mutation profiling, enzyme kinetic assays, and transcriptome analysis of the adapted strain allowed us to decipher the basis of glycerol adaptation at the molecular level. Importantly, increased substrate influx mediated by the mutant glpK and modulation of intracellular cAMP levels were the key drivers of improved fitness in the glycerol-limited condition. Leveraging the enhanced capability of glycerol utilization in the strain, we constructed a GABA-producing E. coli W-derivative with superior GABA production compared to the wild-type. Furthermore, rationally designed inactivation of the non-essential metabolic genes, including ackA, mgsA, and gabT, in the glycerol-adapted strain improved the final GABA titer and specific productivity by 3.9- and 4.3-fold, respectively, compared with the wild-type.

Ageing and rejuvenation models reveal changes in key microbial communities associated with healthy ageing.

BACKGROUND: The gut microbiota is associated with diverse age-related disorders. Several rejuvenation methods, such as probiotic administration and faecal microbiota transplantation, have been applied to alter the gut microbiome and promote healthy ageing. Nevertheless, prolongation of the health span of aged mice by remodelling the gut microbiome remains challenging. RESULTS: Here, we report the changes in gut microbial communities and their functions in mouse models during ageing and three rejuvenation procedures including co-housing, serum-injection and parabiosis. Our results showed that the compositional structure and gene abundance of the intestinal microbiota changed dynamically during the ageing process. Through the three rejuvenation procedures, we observed that the microbial community and intestinal immunity of aged mice were comparable to those of young mice. The results of metagenomic data analysis underscore the importance of the high abundance of Akkermansia and the butyrate biosynthesis pathway in the rejuvenated mouse group. Furthermore, oral administration of Akkermansia sufficiently ameliorated the senescence-related phenotype in the intestinal systems in aged mice and extended the health span, as evidenced by the frailty index and restoration of muscle atrophy. CONCLUSIONS: In conclusion, the changes in key microbial communities and their functions during ageing and three rejuvenation procedures, and the increase in the healthy lifespan of aged mice by oral administration of Akkermansia. Our results provide a rationale for developing therapeutic strategies to achieve healthy active ageing. Video abstract.

Re-classification of Streptomyces venezuelae strains and mining secondary metabolite biosynthetic gene clusters.

Streptomyces species have attracted considerable interest as a reservoir of medically important secondary metabolites, which are even diverse and different between strains. Here, we reassess ten Streptomyces venezuelae strains by presenting the highly resolved classification, using 16S rRNA sequencing, MALDI-TOF MS protein profiling, and whole-genome sequencing. The results revealed that seven of the ten strains were misclassified as S. venezuelae species. Secondary metabolite biosynthetic gene cluster (smBGC) mining and targeted LC-MS/MS based metabolite screening of S. venezuelae and misclassified strains identified in total 59 secondary metabolites production. In addition, a comparison of pyrrolamide-type antibiotic BGCs of four misclassified strains, followed by functional genomics, revealed that athv28 is critical in the synthesis of the anthelvencin precursor, 5-amino-3,4-dihydro-2H-pyrrole-2-carboxylate (ADPC). Our findings illustrate the importance of the accurate classification and better utilization of misclassified Streptomyces strains to discover smBGCs and their secondary metabolite products.

Different Regulatory Modes of Synechocystis sp. PCC 6803 in Response to Photosynthesis Inhibitory Conditions.

Cyanobacteria are promising industrial platforms owing to their ability to produce diverse natural secondary metabolites and nonnative value-added biochemicals from CO2 and light. To fully utilize their industrial potency, it is critical to understand their photosynthetic efficiency under various environmental conditions. In this study, we elucidated the inhibitory mechanisms of photosynthesis under high-light and low-temperature stress conditions in the model cyanobacterium Synechocystis sp. PCC 6803. Under each stress condition, the transcript abundance and translation efficiency were measured using transcriptome sequencing (RNA-seq) and ribosome profiling, and the genome-wide transcription unit architecture was constructed by data integration of transcription start sites and transcript 3'-end positions obtained from differential RNA-seq and sequencing of 3'-ends (Term-seq), respectively. Our results suggested that the mode of photosynthesis inhibition differed between the two stress conditions; high light stress induced photodamage responses, while low temperature stress impaired the translation efficiency of photosynthesis-associated genes. In particular, poor translation of photosystem I resulted from ribosome stalling at the untranslated regions, affecting the overall photosynthetic yield under low temperature stress. Our comprehensive multiomics analysis with transcription unit architecture provides foundational information on photosynthesis for future industrial strain development. IMPORTANCE Cyanobacteria are a compelling biochemical production platform for their ability to propagate using light and atmospheric CO2 via photosynthesis. However, the engineering of strains is hampered by limited understanding of photosynthesis under diverse environmental conditions such as high-light and low-temperature stresses. Herein, we decipher the transcriptomic and translatomic responses of the photosynthetic efficiency to stress conditions using the integrative analysis of multiomic data generated by RNA-seq and ribosome profiling, respectively. Through the generated massive data, along with the guide of the genome-wide transcription unit architecture constructed by transcription start sites and transcript 3'-end positions, we identified the factors affecting photosynthesis at transcription, posttranscription, and translation levels. Importantly, the high-light stress induces photodamage responses, and the low-temperature stress cripples the translation efficiency of photosynthesis-associated genes. The resulting insights provide pivotal information for future cyanobacterial cell factories powered by the engineering toward robust photosynthesis ability.

RiboRid: A low cost, advanced, and ultra-efficient method to remove ribosomal RNA for bacterial transcriptomics.

RNA sequencing techniques have enabled the systematic elucidation of gene expression (RNA-Seq), transcription start sites (differential RNA-Seq), transcript 3' ends (Term-Seq), and post-transcriptional processes (ribosome profiling). The main challenge of transcriptomic studies is to remove ribosomal RNAs (rRNAs), which comprise more than 90% of the total RNA in a cell. Here, we report a low-cost and robust bacterial rRNA depletion method, RiboRid, based on the enzymatic degradation of rRNA by thermostable RNase H. This method implemented experimental considerations to minimize nonspecific degradation of mRNA and is capable of depleting pre-rRNAs that often comprise a large portion of RNA, even after rRNA depletion. We demonstrated the highly efficient removal of rRNA up to a removal efficiency of 99.99% for various transcriptome studies, including RNA-Seq, Term-Seq, and ribosome profiling, with a cost of approximately $10 per sample. This method is expected to be a robust method for large-scale high-throughput bacterial transcriptomic studies.

Genome-Scale Analysis of Acetobacterium woodii Identifies Translational Regulation of Acetogenesis.

Acetogens synthesize acetyl-CoA via the CO2-fixing Wood-Ljungdahl pathway. Despite their ecological and biotechnological importance, their translational regulation of carbon and energy metabolisms remains unclear. Here, we report how carbon and energy metabolisms in the model acetogen Acetobacterium woodii are translationally controlled under different growth conditions. Data integration of genome-scale transcriptomic and translatomic analyses revealed that the acetogenesis genes, including those of the Wood-Ljungdahl pathway and energy metabolism, showed changes in translational efficiency under autotrophic growth conditions. In particular, genes encoding the Wood-Ljungdahl pathway are translated at similar levels to achieve efficient acetogenesis activity under autotrophic growth conditions, whereas genes encoding the carbonyl branch present increased translation levels in comparison to those for the methyl branch under heterotrophic growth conditions. The translation efficiency of genes in the pathways is differentially regulated by 5' untranslated regions and ribosome-binding sequences under different growth conditions. Our findings provide potential strategies to optimize the metabolism of syngas-fermenting acetogenic bacteria for better productivity. IMPORTANCE Acetogens are capable of reducing CO2 to multicarbon compounds (e.g., ethanol or 2,3-butanediol) via the Wood-Ljungdahl pathway. Given that protein synthesis in bacteria is highly energy consuming, acetogens living at the thermodynamic limit of life are inevitably under translation control. Here, we dissect the translational regulation of carbon and energy metabolisms in the model acetogen Acetobacterium woodii under heterotrophic and autotrophic growth conditions. The latter may be experienced when acetogen is used as a cell factory that synthesizes products from CO2 during the gas fermentation process. We found that the methyl and carbonyl branches of the Wood-Ljungdahl pathway are activated at similar translation levels during autotrophic growth. Translation is mainly regulated by the 5'-untranslated-region structure and ribosome-binding-site sequence. This work reveals novel translational regulation for coping with autotrophic growth conditions and provides the systematic data set, including the transcriptome, translatome, and promoter/5'-untranslated-region bioparts.

Multi-Omic Analyses Reveal Habitat Adaptation of Marine Cyanobacterium Synechocystis sp. PCC 7338.

Cyanobacteria are considered as promising microbial cell factories producing a wide array of bio-products. Among them, Synechocystis sp. PCC 7338 has the advantage of growing in seawater, rather than requiring arable land or freshwater. Nonetheless, how this marine cyanobacterium grows under the high salt stress condition remains unknown. Here, we determined its complete genome sequence with the embedded regulatory elements and analyzed the transcriptional changes in response to a high-salt environment. Complete genome sequencing revealed a 3.70 mega base pair genome and three plasmids with a total of 3,589 genes annotated. Differential RNA-seq and Term-seq data aligned to the complete genome provided genome-wide information on genetic regulatory elements, including promoters, ribosome-binding sites, 5'- and 3'-untranslated regions, and terminators. Comparison with freshwater Synechocystis species revealed Synechocystis sp. PCC 7338 genome encodes additional genes, whose functions are related to ion channels to facilitate the adaptation to high salt and high osmotic pressure. Furthermore, a ferric uptake regulator binding motif was found in regulatory regions of various genes including SigF and the genes involved in energy metabolism, suggesting the iron-regulatory network is connected to not only the iron acquisition, but also response to high salt stress and photosynthesis. In addition, the transcriptomics analysis demonstrated a cyclic electron transport through photosystem I was actively used by the strain to satisfy the demand for ATP under high-salt environment. Our comprehensive analyses provide pivotal information to elucidate the genomic functions and regulations in Synechocystis sp. PCC 7338.

Elucidating the Regulatory Elements for Transcription Termination and Posttranscriptional Processing in the Streptomyces clavuligerus Genome.

Identification of transcriptional regulatory elements in the GC-rich Streptomyces genome is essential for the production of novel biochemicals from secondary metabolite biosynthetic gene clusters (smBGCs). Despite many efforts to understand the regulation of transcription initiation in smBGCs, information on the regulation of transcription termination and posttranscriptional processing remains scarce. In this study, we identified the transcriptional regulatory elements in ß-lactam antibiotic-producing Streptomyces clavuligerus ATCC 27064 by determining a total of 1,427 transcript 3'-end positions (TEPs) using the term-seq method. Termination of transcription was governed by three classes of TEPs, of which each displayed unique sequence features. The data integration with transcription start sites and transcriptome data generated 1,648 transcription units (TUs) and 610 transcription unit clusters (TUCs). TU architecture showed that the transcript abundance in TU isoforms of a TUC was potentially affected by the sequence context of their TEPs, suggesting that the regulatory elements of TEPs could control the transcription level in additional layers. We also identified TU features of a xenobiotic response element (XRE) family regulator and DUF397 domain-containing protein, particularly showing the abundance of bidirectional TEPs. Finally, we found that 189 noncoding TUs contained potential cis- and trans-regulatory elements that played a major role in regulating the 5' and 3' UTR. These findings highlight the role of transcriptional regulatory elements in transcription termination and posttranscriptional processing in Streptomyces sp.IMPORTANCE Streptomyces sp. is a great source of bioactive secondary metabolites, including antibiotics, antifungal agents, antiparasitic agents, immunosuppressant compounds, and other drugs. Secondary metabolites are synthesized via multistep conversions of the precursor molecules from primary metabolism, governed by multicomplex enzymes from secondary metabolite biosynthetic gene clusters. As their production is closely related with the growth phase and dynamic cellular status in response to various intra- and extracellular signals, complex regulatory systems tightly control the gene expressions related to secondary metabolism. In this study, we determined genome-wide transcript 3'-end positions and transcription units in the ß-lactam antibiotic producer Streptomyces clavuligerus ATCC 27064 to elucidate the transcriptional regulatory elements in transcription termination and posttranscriptional processing by integration of multiomics data. These unique features, such as transcript 3'-end sequence, potential riboregulators, and potential 3'-untranslated region (UTR) cis-regulatory elements, can be potentially used to design engineering tools that can regulate the transcript abundance of genes for enhancing secondary metabolite production.

Mini review: Enzyme-based DNA synthesis and selective retrieval for data storage.

The market for using and storing digital data is growing, with DNA synthesis emerging as an efficient way to store massive amounts of data. Storing information in DNA mainly consists of two steps: data writing and reading. The writing step requires encoding data in DNA, building one nucleotide at a time as a form of single-stranded DNA (ssDNA). Once the data needs to be read, the target DNA is selectively retrieved and sequenced, which will also be in the form of an ssDNA. Recently, enzyme-based DNA synthesis is emerging as a new method to be a breakthrough on behalf of decades-old chemical synthesis. A few enzymatic methods have been presented for data memory, including the use of terminal deoxynucleotidyl transferase. Besides, enzyme-based amplification or denaturation of the target strand into ssDNA provides selective access to the desired dataset. In this review, we summarize diverse enzymatic methods for either synthesizing ssDNA or retrieving the data-containing DNA.